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ABSTRACT 

Accelerated  life  testing  of  a  product  is  commonly  used  to 
estimate  parameters  of  devices  with  extremely  long  lifetimes.  The 
problem  of  analyzing  accelerated  life  tests  when  the  pioduet  can 
fail  from  any  one  of  p  independent  causes  is  considered.  For  each 
component  it  is  assumed  that  the  time  to  failure  distribution  is 
exponential  with  a  hazard  rate  which  follows  a  "power  rule."  The 
method  of  maximum  likelihood  is  used  to  obtain  estimators  of  i  lie 
power  rule  parameters  when  the  life  test  is  either  type  T,  type  II, 
or,  progressively  censored.  Estimators  of  each  component's  mean 
failure  time  and  survival  function  are  obtained.  Estimators  of 
system  or  subsystem  survival  functions  are  also  obtained.  Results 
of  a  simulation  study  are  presented. 


1  .  rn  t  i  <  nine  t  i  <  <n 

Accelerated  life  lasting  of  a  product  is  ofl.i-n  us.  *d  to 
obtain  i  n  on  mat  i  on  on  its  performance  under  normal  use  roml  i  i  ion.;. 
Ruoh  lasting  involves  sub  joe  1  i  ng  tost  items  to  conditions  more 
extreme  than  usually  encountered  in  normal  use.  This  results  in 
doc reasi ng  the  item’s  mean  life  and  reduces  test  time  and  costs. 
The  estimates  obtained  at  these  severe  conditions  are  extrapolated 
to  no  iina  1  use  conditions.  This  procedure  is  summarized  in  Chapter 
9  of  Mann,  Schafer,  and  R i nypu rwa 1  la  (1  974). 

Little  work  has  been  done  on  analyzing  accelerated  life 
tests  when  an  item  can  fail  from  any  one  of  several,  dist  i  net 
e. sir.es.  Suppose  an  item  consists  of  p  components  such  that  the 
failure  of  any  one  component  causes  the  item  to  fail.  For  each 
ili'in  data  consist  of  the  time  at  which  the  item  fails  and  know¬ 
ledge  of  which  component  failed.  The  problem  is  to  estimate  the 
painmoters  of  the  failure  distribution  of  each  component  at  normal 
use  conditions  from  data  collected  in  an  accelerated  life  test. 
When  !  h>  ic  are  data  at  use  e. unlit  iens  only,  this  problem  Inis  been 
cons  i doped  by  David  and  Moe:..  Lbei  ger-  (197  8  ),  and  Klabf  leisch  and 
:  '  •  nt  i  e  c*  (1980  ).  Del:-, .in  (  1  *J  7  M  )  has  .  .  •i.s  i  doped  <  he  .  u-e.  •  1 .  •  ta  t  ed 
>  i  *  •  test  piiblem  with  comj-'t  i  ng  ri  ;ks  f<-.>r  normal  populations. 

Formally,  let  V^,...,V  be  s  stress  levels  at  which  the 
icci'l  orated  life  test  is  conducted.  At  stress  V.,  let  X.  .  denote 

j.  ^ 

he  1 i fet  i me  of  the  j  component ,  j  =  1 , . . . ,p,  i  =  1 , . . . , s . 

■  ;  ,  |::,I  •  that  t  lie  o  mi  pom  -11 1  life  lengths  aie  i  Iidepe  nden  t  and  each 

lifetime  is  exponent  ini ly  distributed  with  survival  function 

i’  (X.j  ■>  >:)  •  IV  ^  (x)  "Xp(-h  -x),  x  ■  0,  ).  ij  >  0,  (1.1) 

i  -=  1 ,  .  .  .  ,  s 
J  = 


1  »  •  •  .  »P 


Assume  lhat  the  stress  and  component  life  distributions  are  rol  ii  >d 
by  the  power  rule  model. 


Xi j  "  exp (AjJ Vi  J ' 


(1.2) 


where  A.,  B.  axe  constants  to  be  estimated  from  the  accelerated 
3  3 

life  test. 

In  Section  2  the  method  of  maximum  likelihood  is  used  to  ob¬ 
tain  estimates  of  and  for  various  censoring  schemes.  In 
Section  2.1  type  I  censoring  is  considered,  Section  2.2  considers 
type  II  censoring,  while  Section  2.3  discusses  progressive 
censoring . 

The  estimators  of  A^  and  B^  obtained  in  Section  2  are  used 
in  Section  3  to  obtain  estimators  of  component  mean  survival  time 
under  normal  use  conditions.  Confidence  intervals  and  reduced 
bias  estimators  of  these  parameters  are  also  obtained.  Estimates 
of  the  component  survival  function  and  system  survival  function 
are  obtained.  The  problem  of  estimating  system  life  when  a  group 
of  components  are  redesigned  to  increase  reliability  is  also 
considered . 

In  Section  4  results  of  a  simulation  study  are  presented. 

2.  Maximum  Likelihood  Estimation  of  the 
—  Power  ~Ru.Te"  "Parame'Eers' 

Let  V, ,  V_ , ,  V  be  the  stress  levels  at  which  the  accol- 
1  2  s 

crated  life  test  is  to  be  conducted.  At  each  stress  level  we 
assume  that  the  item  survives  until  one  of  its  p  components  fails. 


ho*-  X.^,  denote  the  life  lengths  of  the  p  components 

of  an  item  on  test  at  stress  V. .  Assume  that  the  X. .'s  are 

i  13 


3 


independent  exponential  random  variables  with  hazard  rate  (1.2). 


B  . 

3 


=  vp  J  exp ( A  . ) ,  j  =  l,...,p,  i  =  l,...,s. 


(2.1) 


For  an  item  put  on  tost  at  stress  \A  we  observe 

Y.  =  minimum  (X. .,..., X.  )  and  an  indicator  variable  which  de- 
x  1 1  i  p 

scribes  which  of  the  p  components  failed.  The  method  of  maximum 

likelihood  is  used  to  find  estimators  of  A.  and  B.  for  various 

3  3 

censoring  schemes. 

2 . 1  Type  I  Censoring 

At  each  of  the  s  stress  levels,  n.  items  are  put  on  test, 

l 

i  -  l,...,s.  Testing  is  continued  until  some  fixed  time  . 

Those  r.'s  may  be  different  for  each  stress  level  to  allow  for 

increased  testing  time  at  low  stresses.  Such  a  testing  scheme  is 

a  type  I  censored  accelerated  life  test. 

At  stress  \A  suppose  items  have  failed  prior  to  time 

T  .  ,  i  --  1 ,  .  .  .  ,  s  .  Let  Y.^,  Y  ^  2  '  •  •  •  '  Yb  r  denote  the  corresponding 
1  l 

failure  times.  That  is,  Y^  =  min(X^^^,  x^2k '  ’  *  *  '  Xipk^  w^ere  X^k 
is  the  lifetime  of  the  j*"'"1  component  of  the  k*1^1  item  which  failed 
prior  to  time  t  at  stress  VA  ,  j  =  l,...,p,  k  =  l,...,r^, 
i  =  l,...,s.  Let  m.  .  denote  the  number  of  items  which  failed  due 

13  P 

to  failure  of  component  j,  j  =  l,...,p.  Note  that  r.  -  £  m. . 

1  •  A  J 

3  =  1 

Define  the  total  time  on  test  by 


r  ^ 

Ti  "  kE=1Yik  +  (ni  '  ri)Ti'  1  “  1 - P- 


(2.2) 


The  1  ikel  ihood  function  of  interest  is  i  ■  »  :  i  or. a i  *  o  ‘  he 

•  i  h-  :  ’  •  ih  :  -it  .  *  ;  •  level.  That  is  , 


•5 


s 

l  “  n 

i=l 


p  m .  . 

n  X.  oxp ( -X .  . T . ) 

in  t  i  ■»  ' 

j=l 


tj  i 


(2.3) 


Interchanging  the  order  of  multiplication  we  see  that  L  factors 
into  the  product  of  p  likelihoods,  one  for  each  component.  That 
is , 


L  a  11  I.  whore 

j-1  3 


(2.4) 


s  m.  . 
1 J 


L.  <*  fl  X .  .  J  exp  (-X  .  .T.  )  ,  i  =  l,...,p. 
3  i=1  ID  ID  i  1 


Hence,  it  suffices  to  maximize  each  comj:>onent  likelihood  L.  indi- 

A  A 

vidually.  The  resulting  maximum  likelihood  estimators,  (A.,  B.) 

A  A  " 

will  be  asymptotically  independent  of  (A^,  B^)  when  j  f  kas  all 
The  likelihood  equations  arc 


A .  --In 
1 


(iro)  -  ,n(j1TiviBj) 


j  =  i , . . . ,p. 


(2.5) 


and  I  E  m.  InV. 
a  =  l  ^ 


j  -  1, .  .  .  ,p 


(2.6) 


Bj  can  be  obtained  by  solving  (2.6)  numerically  using  the  Nowton- 
Raphson  procedure.  Clearly  A ^  is  determined  by  (2.5)  once  is 
found . 

Fiora  (2.4)  ,  the  second  derivitives  are  found  to  be 


S  fn  L  .  s  B  . 

2  3  =  -  Z  Tivi  31  oxp(A  )  ,  j  -  l,...,p, 

X A  i  =1  ’ 

1 


(2.7) 


6  Cn  L. 

6  A  .  SB  . 

3  3 


£  T  V.  3  exp (A . ) Cn  V.  ,  j  =  1,... 
i-1  1  1  3  1 


(2.8) 


6  2  f  n  L  . 
3 

B  ,2 
3 


SB. 

>:  T.  V.  3  exp  (A  . )  (InV.)  ,  j  -  1 ,  .  . . ,p.  (2.9) 

i  =1  *  L  3  1 


To  find  the  information  matrix,  note  that  Y^k  has  an  oxpon* 

P 

neniial  distribution  with  hazard  rate  X .  -  Z  X... 

x  j  -i  *3 

Writing,  E(Y;k)  *  RCfJZ^  =--  l)P(Z.k  =  1)  +  Effjjz.^  =  0)P(Z.R  ■---  0), 
wh' ■  i o  Z^k  -  1  if  Y.k  <  t  ^  and  Z^k  =  0  if  Y^k  >  r  ^ ,  it  ran  be  shown 
that  K('Ih)  -  n^n^/X^  where  n ^  —  1  —  exp(-Xar^),  i  -  l,...,s 
is  the  probability  an  item  will  fail  before  time  t  .  The  inverse 
of  the  information  matrix  is 


s  X  .  .n  .  !1 .  ~ 

Z  .J-JLA.Af tnv.)2 
i-1  X.  1 


z .  --  Tr- 

3  °3 


s  X  .  .  n  .  II  . 
-  Z  13  1  L- 
i  -1  X, 


5-n  V. 


s  X .  .  n . n  .  t  nV . 

_  z  J-XJ-.  A  1 

i=l  Xi 

s  X  .  .  n  .  7i . 

Z  13  a  l 
i=l  X; 


win- re 


ss  X  .  .  X,  .  n  .  n,  n .  n. 

D  =  Z  E  --2 - L1--  AJiinv  tn(V./V  ),  j  =  l,...,p.  (2.10) 

3  iA  x.xk  1  1  k 


From  David  and  Mocschbc-rger  ( 1  9  7  8)  note  that  the  im  x  imum 

m .  .  r  . 

likelihood  estimator  of  X. .  is  m  3  and  of  X.  is  .  A  consistant 

l]  T.  l  T. 

•  J  1  l 

:  m.i i  or  of  Z  .  is 
3 


t> 


E  m .  .  ( £n  V .  )  -  I  m.  .  !nV. 

i=l  ^  1  i=i  13 


j  -  E  m .  .  £n  V  . 
J  ..,13  i 


s 

Em.. 
i  =  l  ^ 


whore 


D  -  EEm .  m  tn  V.  ln(V./V  ),  j  = 
3  i^k  13  *3  1  1  k 


(2.11) 


2.2  Type  II  censoring 

For  typo  IT  censor  i  ng  testing  coni  i  lines  until,  a  in  eassigned 
number  of  failures  is  observed.  At  each  stress  level,  V^, 
i  -  l,...,s,  n^  items  are  put  on  tost.  Let  Y^  denote  the  failure 


time  of  the  k  item  put  on  test,  that  is,  Y. 


mi  n  (X.  .,) 
j=l , . . . , p  ^ 


t  h  th 

who  re  x^jk  is  the  life  length  of  the  j  component  of  the  k  item 
put  on  tost  at  stress  VL  ,  j  =  l,...,p,  k  =  l,...,n^,  i  =  l,...,s. 
For  this  testing  scheme,  testing  at  stress  VL  continues  until  a 
preassigned  number,  r^,  of  items  has  failed.  Let  Y^ (p)  <  • • • (r  ^ 

be  the  ordered  failure  times  of  the  r.  fa ;  lures.  Let  m.  .  denote 

i  il 

‘he  number  of  items  which  fail  due  to  failure  of  component 

j,  j  =  l,...,p.  For  type  II  censoring  the  total  time  on  test  is 

b  -  k^Yuk)  *  lni  -  ri>  hep  - 1  -  1 . s-  (2-,2> 

The  likelihood  function  of  interest  for  type  II  censoring  is  iden¬ 
tical  to  (2.3)  with  T^  as  in  (2.11).  Hence,  and  can  be  ob¬ 
tained  as  in  Section  2.1. 

To  find  the  in  formation  natrix,  note  that  the  second  deriva¬ 


tives  of  I,_.  are  as  in  2.7,  2.8,  and  2.9.  l-Y-nce 


1 


s  X  .  . 

"  ^r.tfnvY 

1  i-1  i  1  1 


•  ).  ,  J-  r.  hi  V. 

•  -i  X  .  i  1 1 

l-l  i 


-  I  T J-  r  .  In  V. 
i-lAi  1 


s  X  .  . 
t  t-3 
l  ,  •  -  r  . 

.  ,  X  .  i 

l-l  l 


who  rc 


X  ■  ■  X.  .  r  .  r, 
i]  k]  l  k 


2n  V.  en(V./Vk) ,  j  -  1, 


(2.13) 


The  maximum  likelihood  estimators  of  X.  and  X. .  are 

i  il 

^  m .  .  ~  r . 

X  .  .  '  -1  ^  and  X.  ;  1  ,  respectively .  Substi  Luting  i  he  so  values 
11  'Ih  i  Ti  J 

in  (2.12)  I  ho  maximum  1  ikelihood  estimator  of  Y.  ^  is  obtained. 

This  estimator  is  given  by  (2.11)  . 

2.3  Progressive  Censoring 

Progressive  censoring  is  a  useful  scheme  to  reduce  test 
time  and  cost  while  still  including  some  observations  with  extreme¬ 
ly  long  life  times.  Tn  this  scheme,  at  several  fixed  times  some 
fixed  number  of  items  are  removed  from  study.  This  scheme  is  dis¬ 
cussed  in  Cohen  (1963)  for  the  case  of  normal  and  exponential  data 
from  a  population  with  a  single  failure  mode  when  there  is  not  an 
.iceol  erated  life  lest. 


At  stress  level  V ^ ,  i  -  l,...,s  suppose  that  items  are 
put  on  test.  Censoring  will  occur  in  stages.  At  times 

Til'  1  i  2  '  '  "  ‘  '  T  iM.  -  1  a  f  i  xr-r^  number  r.}.  items  are  removed  from 
the  accelerated  life  test. 

At  time  to  1  he  test  is  either  terminated  with  a  i andom 
i 

tar  :n*r  r.,^  items  still  functioning  or  a  fixed  number  items 

are  t-  moved.  We  assume  that  N ^  is  sufficiently  large  so  that 


b 


there  is  at  least  r  ,  k  =  1,...,M.  items  surviving  to  lime  t  to 

1  K  1  1  K 

Ire  censored.  T.et  m.  .  Ire  the  number  of  items  which  fail  due  i  o 

1 3  p 

failure  of  component  j,  j  =  and,  let  n.  -  >1  m .  .  be  the 

1  j  -- 1  1 3 

total  number  of  failures.  let  Y Y .  be  the  failure  times 

1 1  in. 

i 

of  the  n .  failures.  As  before  ,  Y . ,  -  mi n (X X .  ,  )  whore  X. 

i  lk  ilk  ipk  ljk 

is  the  lifetime  of  the  component  of  the  k*1'"1  item  which  fails 

at  stress  \b  ,  j  =  l,...,p;  k  -  l,...,n^,  i  -  l,.,.,s.  For  this 
scheme  the  total  time  on  test  is  defined  by 


n  .  m  . 

l  I 

T .  =  EY..+  E  r  . ,  r  • ,  ,  i  -  1 , 

1  k=l  lk  k-1  ,k  lk 


,  s . 


(2.11) 


As  before  it  can  bo  shown  that  the  likelihood  of  interest  is  as  in 

(2.4).  The  maximum  likelihood  estimators  of  A.  and  B.  as  well  as 

11 

the  second  derivatives  are  as  in  Section  2.1. 

For  this  censoring  scheme  the  information  matrix  and  an 
estimator  of  the  matrix  are  obtained  in  the  Appendix. 

3.  Estimation  of  population  parameters  at  use  conditions 


Suppose  an  accelerated  life  tost  has  been  conducted  by  one 


of  the  censoring  schemes  discussed  in  Section  2.  Let  A.  and  B. 

1  1 

be  the  maximum  likelihood  estimators  of  A.  and  B.,  j  =  l,...,p. 

1  1  J  1 

T.et 


be  Hie  iiiv..  i  'if;  of  informal  on  mat 
A.  10.  T.et  be  the  correspond  i  na 


r  i  x  o i  > '  a  i  n oc  f  i  orr,  (2.10),  (  2  .  )  3  )  nr 

estimator  of  I.  obta i nod  from 


(2.11),  A. 11 ,  or  A. 12. 


Rased  on  this  in  fun.i.it  ion  we  wish  to 


estimate  each  component's  mean  survival  lime  and  suivival  function 

at  the  use  conditions.  Tn  addition  we  shall  obtain  estimators  of 

the  survival  function  for  I  he  entire  sysl.m  of  components  or  a 

subsystem  of  components  under  use  conditions. 

I>et  V  denote  the  s(  i  ess  under  use  condi  1  ions.  bet  X  .  and 
u  u  j 

Uu  j  denote  the  hazard  rate  and  moan  suivival  I  i  me  at  this  si  i.e.s. 
By  (2.1) 


Xuj  =  Vu  ]exPAjf  j  -  1 . p 


(3.2) 


%  ^  Vu  ]exp(-Ad,  j  l - P- 


(3.3) 


By  the  invariance  principle  of  maximum  likelihood  estimators  the 

maximum  likelihood  estimators  of  X  .  and  u  .  are 

uj  uj 


X  B . 

uj  =  vuj  exp (A j ) ,  j  =  l,...,p. 


(3.4) 


u  u j  =  Vu  exp ( -A ^ ) ,  j  =  l,...,p. 


(3.5) 


respectively.  For  large  sample  sizes  at  all  the  stress 

levels  ,  (A^,  )  are  asymptot ical ly  bivariate  norm, ally 

distributed  with  mean  vector  (A.,  B.)  and  variance  l . ,  j  =  l,...,p. 

3  3  3  *  ^ 

y\  A  a 

Since  In  X  .  is  linear  in  A.  and  B.,  it  has  an  asymptotic  no;  ml 
u3  3  3 


:  ’  1  a  :  i  i:  v: :  *  h  r.  in  £  n 


V  r;  f  .  .  1 1  \  CC 


hi’  *  2’i23>  '»V 


(3.61 


1  u 

A 

Similarly,  £n  pujis  asymptot ically  normal  with  mean  £npuj  and 
„  2 

variance  o  .  . 

uj 

By  using  the  properties  of  the  log  normal  distribution  one 
can  show  that  asymptot i eally 


E<Muj)  =  %j  GXP(°uj/2)»  j 


(3.  7) 


and  the  mean  squared  error  of  ij^.  is  OKu^  -  mu^)?)  is 

m.s.e.  (p  .)  =  u  f.  (exp  (2o  2.)  -  2  exp(of./2)  1],  j  =  l,...,p,  <3. 

uj  uj  uj  uj 


with  similar  expressions  for  A 


unbiased  estimator  of  p  .is 

u  J 


ui 


From  (3.7)  we  see  that  an 


p  .  p  .  exp  (-o  ,/2),  j  - 

uj  uj  1  u]  '  J  r 


(3.9) 


which  has  mean  squared  error 


m.s.e.  (Puj)  =  p^(exp(o^)  -  1),  j  =  1, 


(3.10) 


which  is  always  smaller  than  the  mean  squared  error  of  p^  .  Since  e* 

is  unknown  a  reduced  bias  estimater  of  p  .  is  obtained  by  substitul 
«  u  J 

in  (3.9).  Similar  results  hold  for  estimating  A^-. 

The  asymptotic  normality  of  £n  .  c<\n  be  used  to  construct 

(1  -  a)  x  100%  confidence  intervals  for  p  ..  These  are 

uj 


(p  .  exp(-Z,  .),  p  .  exp(Z,  ,~o  .)  ,  j  =  l,...,p 

Muj  ^  l-a/2  uj  uj  1-aA  uj  J  r 

th 


(3.11) 


where  z^_a is  the  1  -  a/2  percentage  [joint  of  a  standard  nm  iml 

random  variable.  A  similar  interval  can  be  obtained  for  >  .. 

u  1 

Let  F  . (t)  -  exp(-tA  .)  be  the  survival  function  of  comoenont 
uj  '  r  uj 

j ,  j  -  l,...,p.  The  maximum  likelihood  estimator  of  F^ft)  is 


Fuj(t)  -  «P(-tXujl,  t  >  0,  j  - 


(3.12) 


2 


This  is  a  conservative  interval  in  the  following  sense.  From  (3.1 i). 


1/P  //  exp  ( Zfto  )  _  _  exp  ( 

(1  -  a)  =  P  F  .  (t)  p  U]  <  F  .  (t)  <  F.  .  ... 

\  uj  -  uj  13 (t) 


V-i’).  fur  j  ’  1  , 


Since  (Fhj(t),  Fu^  (t) )  are  asymptotically  independent  for  j  /  j', 

/-  exp{z  o  .)  _  _  exp(-  Z  o  ) 

1  -  a  -  p(Fuj(t)  3  <  Fuj(t)  <  Fuj(t)  3  ;  for  all  j  cK 


<  p  n  F  •  (t) 

~  WiK  U3 


exp  (Zc°  • )  _  (K)  ~  c: 

t)  3  £  f  (t)  <  n  f  .  (t) 

u  “  jeK  U3 


exp ( -  ZB0uj) 


4.  Simulation  Study 

A  simulation  study  was  conducted  to  compare  the  effects  of 
censoring  on  the  estimation  procedure.  A  two  component  series 
system  of  independent  exponential  component s  was  used.  An  acceler¬ 
ated  life  test  was  conducted  at  3  stress  levels,  =  10,  =  15, 

-20.  At  each  stress  level  the  observations  were  generated  from 
an  exponential  distribution  with  hazard  rate  A^,  i  =  1,...»3 
j  -  1,2,  where  A  ^  is  given  by  (1.2)  with  A^  =  -1,  =  1,  A2  =  -3, 

and  B2  =  2  throughout  the  study.  The  values  of  A_,  and  Hj  =  P 
(component  j  fails  first)  are  given  in  Table  I  for  V^,  and  a 

usage  stress  of  5. 

Table  1 


V  . 

1 

Xil 

Xi2 

n2 

5 

1.8394 

1.2447 

.  5964 

.  4036 

10 

3.6788 

4  .9787 

.4249 

.5751 

1  5 

5.5182 

11 .2021 

.  3300 

.6700 

20 

7  .  3576 

19.9148 

.2698 

.  7302 

■'or  each  censoring  scheme  500  repetitions  of  the  accelerated 


life  test  were  per To l red .  Maximum  livelihood  estimates  of  A ^ ,  and 


1 


Bj>  j  -  1 1  2  were  found  by  using  the  Mi ‘wt on-Haphson  procedure 
discussed  previously.  The  estimated  covariance  matrix  was  found 
by  appropriate  expression.  Both  maximum  likelihood  and  reduced 
bias  estimators  of  the  component  hazard  rate  and  mean  survival  at 
a  use  stress  of  5  were  calculated.  Tables  2-7  report  for  each 
estimator  its  estimated  moan  (standard  error),  estimated  bias 
( | Bias | /Standard  error),  estimated  variance,  and  the  P -value  of  a 
Kol mogorov-Smi rnov  test  for  normality  for  each  component.  Die 
following  censoring  schemes  were  used: 


Scheme  I  -  complete  failure  times  on  all  items,  sample  sizes 
ni  =  20,  i  =  1,  2,  3. 

Scheme  II  -  complete  failure  times  on  all  items,  sample  sizes 


Scheme  III 
Scheme  IV 
Scheme  V 


Scheme  VI 


=  10,  i  =  1,  2,  3. 

type  II  censoring  =  40,  =  20,  i  =  1,  2,  3. 

type  II  censoring,  n^  =  40,  =  10,  i  =  1,  2,  3. 

type  I  censoring  at  the  same  fixed  point  =  40, 
i  =  1,  2,  3,  r1  =  .08,  x2  =  .04,  x3  *  .025, 

E(ri).  =  19.989,  E(r2)  =  19.507,  E(r3)  =  19.772. 
type  I  censoring  at  the  same  fixed  point  n^  -  40, 
i  =  1,  2,  3,  r1  =  .03,  T  =  .017,  t  3  =  .010, 


E(rx)  =  9.15,  E(r2)  =  9.89,  E(r3)  -  9.55. 

Several  conclusions  may  be  drawn  from  the  study.  First  the 
asymptotic  normality  property  of  the  m.l.e.  estimators  up:  s  i  o 

id  fur  the  is-:  ate  :•  1  e  .  S-c  nnlv,  "die  >e'iu,.vd 


(’St  i  mators  X  .  and 
uj 


perform  as  expected  when  the  estimated 


variance  was  used.  Finally,  the  asymptotic  unbiasedness  of  t  lie 
maximum  likelihood  estimators  seems  to  be  reasons)- le  for  these 
moderate  sample  sizes. 

5.  Dependent  Risks 


The  assumption  of  independence  of  the  times  to  failure  of 
each  component  made  in  the  previous  sections  can  he  replaced  by 
the  assumption  that  the  failure  times  of  the  item  follow  a 
Marshall  and  Olkin  (1967)  multivariate  exponential  distribution. 
For  this  model  simultaneous  failures  of  several  components  is 
possible.  Here  there  are  2^-1  possible  failure  modes  which 
represent  system  failure  from  the  simultaneous  failures  of  any 
subset  of  components.  Each  of  these  failure  modes  is  assumed  to 
have  a  time  to  occurrence  distribution  piven  by  1.1  with  hazard 
rate  1.2.  At  use  conditions  the  survival  function  of  the  time  to 
failure  distribution  of  a  particular  component  can  be  estimated 
as  in  the  previous  section. 
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Appendix.  Calculation  of  the  information  matrix 
fo_r  progressively  censored  data 

The  distribution  of  the  total  time  on  test  given  by 

(2.13)  depends  only  upon  the  distribution  of  the  minimum  of  p 

lifetimes  of  the  components  which  is  exponential  with  parameter 

A  ^ .  To  simplify  the  notation  we  shall  suppress  the  subscript 

i  which  refers  to  the  stress  level  under  consideration.  The 

problem  is  to  calculate  the  expected  value  of  the  total  time  on 

test  T  given  by 

n  M 

E  V  +  Z  rit..  (A.l) 

i=l  i=l 

Let  f^  be  the  number  of  items  which  fail  in  the  interval 

[t.  i  =  1 ,  .  .  .  ,M  where  -  0.  Let  S.  be  the  sum  of  the 

i-l  l  0  l 

f^  failure  times  of  those  items  failing  in  t^)  ,  i  =  1,....M* 

Let  =  1  -  exp(-At^)  and  F\  =  1  -  F^f  i  =  1 ,  .  .  .  ,M .  We  consider 
two  cases. 

Case  1.  rM  fixed 

.  .  M 

Suppose  r  is  fixed  so  that  f ,, , ,  =  N  -  Z  (r.  +  f.)  items 

L  M  M+l  l  l 

fail  in  the  interval  (t  ,  00 )  .  Cohen  (1963)  shows  that 

m 

IN  F1 

i-I  rj 

(N  -  Z  — ) (F.  -  F.  ) ,  i  = 

3=1  F3 

where  F,1  +  ^  =  1. 


2  ,  .  .  .  ,M  +  1 


1  7 


To  find  E(S^)  let  Y^,...,Y^  be  the  failure  Limes  of  the 

i 

items  which  fail  in  the  interval  t^),  so  11iat 


S  . 
l 


f . 

l 

I  Y 
K=1 


iK 


for  i 


1 #  •  •  • /M i 


i 


f  . 

l 


E(S,  )  =  EE  (  I  Y.„J  f,  ) 


K-l 


rK'  iJ 


“  E<fiE(YulTi-l  i  Yil  *  V‘ 


(A. 3) 


Now,  given  that  Y ^ ^  fails  in  [ t  ^  ^ ,  t  ^ )  the  conditional  expectation 
is 


E(YillTi-l  <-  Yil  <  V  “ 


T  . 
1 


yXe  ^^dy/(F^  -  F^_^) 


Ti-1 


(A. 4) 


=  (XTi~lF\-l  ~  ATiFi  +  Fi  '  Fi-l>/^Fi  '  Fi-1> 


Combining  A. 2  and  A. 4  by  A. 3  we  see  that 


E(Si)  - 


N( 


t i F i )  for  i  =  1 


(A. 5) 


i-1  r3 


F.  -  F.  , 

(N  -  Z  )(x.  ,F._,  -  T.F.  +  -1 - =—=■)  for  i  =  +  1 

j=l  F.  1-1  1-111  X 

J  J 


So 


M 


m 


7-(N  -  I  r  -  1  r  T  ) 

A  i-1  1  i-1  1 


M+l 

E  (  Z  S.) 
i-1  1 


(A. 6) 


and 


M 

N  -  E  r  i 

E  (T)  =  —  -1"-1  —  . 

X 


(A. 7) 


Case  2.  random 

As  in  Case  1  the  E(f^)  is  given  by  A. 2  for  i  -  and 

E(S^)  is  given  by  A. 5  for  i  =  1, , . . ,M.  It  remains  to  find  E ( r^) . 
Note  that 


M 


M  M-l 

N  -  E  f .  -  £  r .  . 


i=l 


i=X 


Hence 


M  M-l 

E(rM)  =  N  -  E  E(f  .  )  -  £  r.  , 

"  i=l  1  i=l  1 


which  is,  from  A. 2, 


M-l  r . 

E  ( r  )  =  (N  -  E  - -)  F  . 

M  i=l  K 


(A. 8) 


Combining  A. 2,  A. 5,  and  A. 8  we  obtain 


m/  "Xtm\  m_1 

Em-S(l-e  M)  -  .1^,11  -e 


->(tm  -  Ti>' 


(A. 9) 


The  information  matrix  of  interest  for  the  accelerated  life 


test,  using  the  notation  of  Section  2,  is 


z . 

3 


s  m .  .  ~ 


D  . 
3 


Z  E(n.  )  Un  V.  ) 


i=lni 


s  ra .  .  ~ 

-  Z  E  (n. )  2-n  V. 

i=l  ni  1  1 


s 
-  Z 
i=l 


m .  . 
n . 


l 


E(n.)lnV. 
l  i 


s  rn .  .  /v 
z  1J-  E  (n  .  ) 

.  ,  n .  i 

i=l  l 


whore 


/s  in .  .  rn,.  .  ^  a 

D.  =  Z  Z  --il-iil  E(n,)E(nJtnV.  £n(V  /V  ),  j  =  1,... 
3  i^K  ninK  i  K  i  i  K 


and 


E(n.) 


Mi-1 


N.(l  -  exp  {-  X.xM_))  -  ^ul  1  -  e*P(“  Xi(TiMi 


i 


Table  2 


Censoring  Scheme  I 


Comt-onent  I 


|  T.ne 
Val  ue 
of 

Ta  ra¬ 


ft  -l  -1.0207 

(.07  51) 

P  1  2. ('047 

-  ~  (.365) 

V  (  A  )  4.2543  4. 5357 

,  ,  (.0556) 

V  (  E  )  .6143  .6599 

-  „  ..  (.0085) 

V(  B,  A)  -1  -6073  -1.724 

-  (.02)5) 

C  HA  .6094  .5962 

-U  (.0373) 

A  1.8394  2.5703 

(.1157) 

A  1.8394  1.8275 

-U  (.0827! 

0  .5437  .7902 

(.0410) 

Z  .5437  .5366 

(.0257) 

(3  .6717  .1158 

U  (.0085) 


Ti  ue 
Val  ue 
of 

P.ira- 


.  020  7 

(.22) 
.004  7 
(.13) 

.  301  3 
(5.42) 
.0456 

(5. 31) 
.  1  166 
(5.42) 
.01  32 

(.35) 
.  7309 

(6.32) 
.0119 

(  .  14) 
.  2466 
(6.02) 
.0071 
(.28) 
.0440 
(5.18) 


4.32  3 


1 . 5499 


.0359 


.2324 


.6959 


6.6902 


3.414 


.8392 


.0362 


Component  2 


.0001  2.3750 


.0001  .3232 


,0001  -.8714 


.672  .2189 


.  0001  1.2447 


.0001  1  .2447 


.0001  .8034 


.0001  .8034 


.0001  .3891 


-3.0762 
(.0726) 
2.02/7 
(.0268) 
2.4445 
(.0170) 
.  3322 
(.0022) 
-.8964 
(.0022) 
.1873 
(.03) 
1.511 
(.0514) 
1.2389 
(.0432) 
1 .0363 
(.0339) 
.8293 
(.026) 
.4197 
(.0032) 


.0762  I 

(.98) 

.0227 

(1.04) 

.0695 

(4.95) 

.  0090 
(4.20) 
.0249 
(4.15) 
.0315 
(1.05) 

.  266  3 
(5.18) 
.0058 
(.1336) 
.2329 
(6.88) 
.0259 
(.995) 
.0306 
(9.59) 


2.6367 


.  3  586 


.4  500 


1.3211 


.9315 


.  5737 


.  3385 


.0001 


Table  3 


Censoring  Scheme  II 


Component  1  CY.i  p<  >:iont  2 


'!'r  ue 

True 

Value 

Y.i  1  ue 

of 

of 

Moan 

B1  as 

Var 

TV*  i 

F  l  as 

\  ,ir 

Pci  Tel- 

Pr'i  r  r,- 

r,e  1  or 

.  .  . 

j:.  L*‘I 

■ 

1 

A 

-1.0 

-1.132 

.  1  32 

9.6796 

.968 

-3 

-3.0195 

.0195 

4.8315 

•  »  ’8 

(.1  393) 

(.95) 

(.0983) 

(.2) 

1 

B 

1.0 

1 .0376 

.0  376 

1.4061 

.  880 

2.0 

2.0179 

.0179 

.  1  886 

.9  71  | 

-A 

(.0631) 

(.71) 

(  .0363) 

(.49) 

1 

v(  A) 

8.5087 

10.169 

1.6603 

26.2856 

.0001 

4.7  50 

4.9798 

.  2298 

1.1670 

.0001  j 

/s  ^ 

( . 2293) 

(7.24) 

(.0483) 

(4.76) 

] 

V (  B) 

1.2287 

1.4701 

.2414 

.  5195 

.0001 

.6464 

.6774 

.0310 

.0188 

.0001 

^  A  A 

(.0322) 

(7.27) 

(  .0061) 

(5.08) 

v  (  B ,  A ) 

-3.2150 

-3.8444 

-.6294 

3.6331 

.0001 

L1 . 74  30 

Ll .8270 

.  0840 

.1470 

.0001  I 

( .0852) 

(7.39) 

(.0171) 

(4.9) 

i 

.6094 

.  5379 

.0715 

1.5162 

1  .0 

.  2189 

.  2281 

.0092 

.8198 

.796 

( .0051) 

(1.30) 

(  .0409) 

(  .22) 

1 

X 

1.8394 

3.6161 

1 . 7767 

4  3. 6677 

.0001 

1  .  2447 

1.8955 

.6508 

5. 2024 

.  nool 

u 

(.2959) 

(6.01) 

(.1020) 

(6.38) 

X 

1.8394 

1.7695 

.  0699 

9.8243 

.0001 

1  .  2447 

1  .  2822 

.0375 

2.4498 

.  01)0.1 

u 

( . 1403) 

(.5) 

(.0700) 

(-54) 

i 

u 

.  5437 

3.1513 

2.6076 

2004.65 

.0001 

.  8034 

1.2741 

.4707 

5.8434 

.0001  ! 

U 

(2.0043) 

(1.3) 

(.1081) 

(4.35) 

1 

U 

.  54  37 

.4647 

.0790 

.  3  368 

.  0001 

.8034 

.7718 

.0316 

1  .  1955 

.0001 

'2 

( .0259) 

(3.05) 

(.0489) 

(.65) 

0 

1 . 3435 

1 .6023 

.2588 

.07748 

.0001 

.  7780 

.8536 

.0756 

.  °4  1  3 

.  909  1 

u 

( .03937) 

(6.52) 

(  .0091) 

(8.31) 

j 

.  J 

--  -  -  -J 

. 

.  .  -  - 

) 

Table  4 


Censoring  Scheme  III 


Coniocm  1 


Component  2 


...  -  ■ 

True 

'  . 

True 

' 

Value 

Value 

Nor- 

Kc.r- 

of 

of 

Moan 

Bi  as 

Var 

mal 

Mean 

Bias 

V.ir 

1 

Para- 

Para- 

Pie  ter 

’•v  vei 

'  A 

-1.0 

-  .9954 

.0046 

4.8872 

.97 

.3  . j 

-3.105? 

.  1052 

2.1859 

1.0  ! 

( .0909) 

(  .05) 

(  .0661) 

(1.59) 

B 

1.0 

.9962 

.0  338 

.  7107 

.  38 

2 

2.0236 

.0236 

.2991 

1  .0 

r  - 

( .0377) 

(  .90) 

( .0245) 

(  -97) 

V  (A  ) 

4.25.14 

4.6473 

.  3929 

1.8201 

.0001 

2.375 

2.421 

.0457 

.  1  279 

.0002 

A  - 

(  .0603) 

(6.55) 

(  .0160) 

(2.86) 

!  V  (  B  ) 

.6143 

.6742 

.0599 

.0401 

.0011 

.  3232 

.  3289 

.0057 

.0020 

.001 

A  ^ 

( .0090) 

(6.68) 

( .002) 

(2.85) 

V  (  B  ,  A  ) 

-1.607 

-1.7599 

.1526 

.2661 

.0001 

-.8714 

-  .8877 

.0161 

.0159 

.0016 

A 

( .0231) 

(6.63) 

( .0056) 

(2.86) 

£n  A 

.6094 

.5597 

.0497 

.  7641 

.73 

.2189 

.1517 

.0672 

.  3721 

1  .0 

A 

( .0390) 

(1.28) 

.0273 

(2.49) 

A 

1 .8394 

2.5647 

.  7253 

7.5631 

.0001 

1.244  7 

1 . 3966 

.  1519 

.  7971 

.0001 

(.1230) 

(8.39) 

( .0399) 

(3.806) 

x 

1 .8394 

1.8173 

.0221 

3.7407 

.0001 

1.2447 

1  . 1446 

.1000 

.5594 

.  0001 

A 

(  .0  865 

(  .26) 

(.0335) 

(2.859) 

V 

.  54  37 

.8469 

.  70  32 

.9896 

.0001 

.8034 

1 .0351 

.2317 

.4985 

.0001 

(  .0146) 

(( .82) 

( .0306) 

(7.57) 

u 

.  5437 

.  6694 

.0157 

.2867 

.0001 

.8034 

.8320 

.0286 

.2790 

.  Ouol 

( .02  19) 

(  .66) 

.0236) 

(1.21) 

0 

.6717 

.  7287 

.0570 

.0456 

.0001 

.  3891 

.1458 

.0266 

.0046 

.0001 

(  .0090 

(5.95) 

( .0038) 

(8.87) 

_ _  .. 

_ _ 

_ _ . 

.  _ 

.  .  . 

! 

•1 


Table  & 

Censoring  Sdn-rv1  IV 


Component  1  Covpoiient  2 


v 

True 

Value 

of 

Prii:a- 

ir-ter 

Mean 

8;  as 

Var 

... 

Nor¬ 

mal 

True 

Value 

of 

Rara- 

tv/»  t  c-  r 

Mean 

■  • 

R  j  as 

Va  r 

h - 

| 

:r  -> ) 

A 

-1.0 

-1.3776 

.  3776 

10.4405 

.59 

-3.0 

-2.976 

.024 

4.9545 

1.0 

B 

1.0 

(.1445) 

1.0822 

(2.61) 

.0822 

1.4731 

1.0 

2.0 

.0995 

1.9608 

(..24) 

.0392 

.6  78 

.84 

V  (  A) 

8. 5087 

(.0543) 

11.2201 

(1.51) 

2.7114 

109.319 

.0001 

4.750 

(.0368) 

4.9610 

(1.06) 

.211 

1  .  348 

1 

.0001 

V  (  B  ) 

1 . 2287 

(.4670) 

1.6055 

(5.81) 

.  3768 

1.8410 

.0001 

.6464 

(.0519) 

.6752 

(4.07) 

.0288 

.0211 

1 

.0001 

^  vs  vs 

V(B,  A) 

-3.2150 

(.0607) 

-4.2206 

(6.21) 

1.0059 

14.0076 

.0001 

-1.743 

(.0065) 

-1.8206 

(4.43) 

.0776 

.1675 

.0001 

A 

CnA 

.6094 

(.  1673) 

.  3641 

(6.01) 

.2453 

1.718 

.2607 

.  2189 

(.0183) 

.1797 

(4.24) 

.0392 

.8416 

i 

.909  ! 

U 

A 

1.8394 

(.2653) 

3.0967 

(4.34) 

1.2570 

35. 1877 

.0001 

1.2447 

(.0410) 

1.8137 

(.96) 

.  5690 

5.437 

i 

.0001 

A 

1. 8394 

(.2653) 
1 . 4905 

(4. 34) 

.  3488 

7.2648 

.0001 

1.2447 

(.1040) 

1.2413 

(5.47) 

.0037 

2.873 

.0001  ! 

u 

u 

.  5437 

( .  1205) 
4.7719 

(2.89) 

4.2282 

4159.18 

.0001 

.  8034 

(.0758) 

1.2994 

(.05) 

.  4960 

2.464 

1 

.0001  i 

u 

V 

.  5437 

(2.8840) 

.  5566 

(1.46) 

.0129 

.7334 

.0001 

.8034 

(.0702) 

.8017 

(.71) 

.0013 

.  7480 

1 

>01  1 

u 

'2 

0 

1  .  34  35 

(.0  383) 
1.7932 

(.34) 

.0497 

3.7343 

.0001 

.  7780 

(.0387) 

.8500 

(.03) 

.0720 

.0407 

1 

.<•001  . 

u 

(.0864) 

(.58) 

(.0099) 

(7.27) 

( 

1 

( 

Table  6 


Censoring  Scheme  V 


Component  1  Co:i.;oii"nt  2 


True 

" 

True 

| 

Value 

of 

Ik?  an 

Bi  as 

V^r 

NY>1— 

cal 

Value 

°f 

C  -  an 

Bi  as 

VYi  r 

Nor¬ 
mal  j 

Pci 

ia  .-ter 

. 

Ba  ca¬ 
l'  .ter 

-  .  .  . 

1 

I 

( 

* 

!  A 

-1.0 

-1.1598 

.  1  598 

4.5657 

.862 

-3 

-3.0257 

.0257 

2.7150 

455  1 

B 

1.0 

(.0956) 

1.0425 

U  .67) 
.0425 

.6496 

.947 

2 

(.07  37) 
2.0064 

(.35) 

.0064 

.3713 

<  > 

» 

4. 1700 

(.0360) 

4.7316 

(1.18) 

.5616 

2.1813 

.0001 

2.38  56 

(.0273) 

2.4813 

(.23) 

.0957 

.0233 

.0001  i 

!  V(  B) 

.6031 

(.0662) 

.6830 

(8.47) 

.0799 

.0455 

.0001 

.  3249 

(.0233) 
.  3377 

(4.11) 

.  0128 

.0045 

! 

.0001  ; 

|V(B,  A) 

-1.5765 

(.0095) 

-1.7873 

(8.41) 

.2108 

.  3309 

.0001 

-.8757 

(.3030) 

-.9104 

(4.26) 

.0347 

.0348 

1 

.0001  ! 

r  ha 

.6094 

(.0249) 

.5181 

(8.46) 
.091  3 

.7411 

.779 

.2189 

(.0083) 
.  2036 

(4.18) 

.0153 

.4622 

i 

.5846  ! 

U 

A 

X 

1.8394 

(.0385) 
2. 3313 

(2. 37) 
.4919 

3.2589 

.0001 

1.2447 

(.0304) 

1.5240 

(.50) 

.2793 

1.1280 

i 

.oooi  j 

U 

A 

1.8394 

(.0867) 

1.6728 

(5.67) 

.1666 

2.1204 

.0001 

1.2447 

(.0475] 

1.2561 

(5.88) 
.01  14 

.8315 

.0001 

u 

A 

v, 

.5437 

(.0651) 

.9100 

(2.56) 

.  366  3 

1.514 

.0001 

.80  34 

(.0407] 

.9100 

(.28) 

.  1066 

1.5137 

.Oil'll 

U 

Vl 

.54  37 

(.0550) 

.5807 

(6.66) 

.0370 

.4070 

.0001 

.8034 

(.0550] 

.8231 

(1.93) 

.0197 

.4377 

.  lion) 

u 

* 

0 

.6740 

(.0285) 

.7478 

(1-30) 

.0738 

.0592 

.0001 

.4084 

(.0296] 

.4254 

(.67) 

.017 

.0094 

.0001 

u 

(.0109) 

(6.77) 

(.0043! 

-J 

(3.95) 

.  . i 

Table  7 


2  (> 


Censoring  Scheme  VI 


Ci  ■  ut  1  Component  2 


i  '  ue 

"i’i  ue 

" 

’l 

",.1  u-  ‘ 

Value 

) 

Nor  — 

of 

.or-  | 

or 

Vi  .in 

v  <r 

a  1 

Near. 

Bias 

Var 

n  ? )  j 

l'.u  i- 

Pci  r ci  “ 

| 

■ ..  :  .■  r 

.  . 

uriir 

. 

•  - - 

.....  .  . 

, 

A 

-1  .0 

-1 . 3115 

.3115 

12.522 

.118 

-3 

-3.1109 

.1109 

6.1581 

.1760 

(  .  158  3) 

(1  .97) 

(.1110) 

(1.00) 

! 

a 

1  .0 

1.0715 

.0715 

1.7535 

.218 

2 

2.0223 

.0223 

.8292 

.6622 j 

A 

(  .0592) 

(1.21) 

(.0407) 

(.55) 

V(  \  ) 

8.9  299 

1  3.  3662 

4 . 4  36  3 

181.61? 

.0001 

5.1163 

5.9371 

.8208 

6.8142 

.0001 

)a 

(  .6027) 

(7. 36) 

(.1167) 

(7.03) 

1 

j 

V(  B  ) 

l  .  2809 

1.8959 

.6150 

2.9167 

.0001 

.69  39 

.8018 

.1079 

.1079 

.0001  I 

A 

<  .0764) 

(8.05) 

(.0147) 

(7.34) 

i 

V  ,  fc  ,  A  ) 

-3. 3588 

-5.0066 

-1.6478 

22.7431 

.0001 

-1.8744 

-2.1709 

.2965 

.8511 

.0001 j 

A 

(  .2133) 

(7.73) 

.2189 

( .0413) 

(7.18) 

‘n 

.  8094 

.4130 

.  1964 

2.0741 

.025 

.  1439 

.0750 

1  .0688 

.010 

•- 

.0644 1(3.05) 

( .0422) 

(1.78) 

] 

1  \ 

1  .8393 

3.4528 

1.6135 

41  .  7061 

.0001 

1.24/7 

1 . 7861 

.5414 

2.8247 

.000) 

u 

(  .  20H8) 

(5.59) 

( .0752) 

(7.20) 

>■ 

1 .839  3 

1.5791 

.  2602 

6.9537 

.0001 

1 .2477 

1.4876 

.0571 

1.4911 

.0001 

Au 

(.11  79) 

(2.21) 

( .0546) 

(1.05) 

1 

.-437!  5.4206 

4.8769 

1580 . 64 

.0001 

.8034 

1.7481 

.9447 

19.2456 

.0001 

(  1  .  -’78) 

(2.74) 

( .1962) 

(4.81) 

1 

1  C 

.  5437 

.4494 

.0943 

.  6336 

.0001 

.8034 

.8089 

.0055 

1  .  1809 

.0001 , 

1 

u 

'  A  A 

(  .0356) 

(2.65) 

( .0486) 

(.11) 

! 

1,4362 

2.1615) 

.7253 

6.6171 

.0001 

.  8804 

1.0262 

.  1458 

.2453 

.0001 

1 

1 

(.1150) 

(6.31) 

( .0221) 

(6.63) 

- . -  -- 

1  -  ■  -  '  -  — 

* 

-  - 

.  . 

